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A measurement of the W boson mass using 4.3 fb 1 of integrated luminosity collected with 
the DO detector during Run II of the Fermilab Tevatron collider is presented. Based on 



, a sample of 1.68 x 10 B W — > ev candidate events, the W boson mass is measured to be 

Mw = 80.367 ± 0.026 GeV. Combining this result with an earlier DO result determined 
using an independent Run II data sample of 1 fb 1 of integrated luminosity, yields Mw = 
1.375 ± 0.023 GeV. 



1 Introduction 

The standard model (SM) of particle physics predicts the existence of a hypothetical scalar par- 
ticle, the Higgs boson, as a result of the spontaneous electroweak symmetry breaking mechanism 
that explains the masses of the vector bosons (W and Z). Direct searches ^ of the SM Higgs 
boson have limited its possible mass ranges to be 115-127 GeV or above 600 GeV at 95% C.L.. 

Also predicted by the SM, is a relationship between the W boson mass (M^) and the Higgs 
boson mass (Mh), together with other electroweak parameters such as the top quark mass (iWt). 
Precise knowledge of the value of My/ (and other parameters), therefore, can be used to predict 
the possible mass range of the hypothetical Higgs boson. 

A combination^ of previous measurements of My/ before winter 2012 yielded a world average 
value of My/ = 80.399 ± 0.023 GeV. This result together with the current Mt measurement^ 
(and measurements of other electroweak parameters) predicts^ a mass range of the Higgs boson 
of M H = GeV, with an upper limit of 161 GeV at 95% C.L.. The predicted M H range 

overlaps with the ranges allowed by direct searches. However, the predicted range is large, and 
the experimental precision on the My/ is so far the limiting factor in this prediction. Therefore, 
improving the precision of the My/ measurements to narrow the predicted Mh range, and 
comparing the prediction with direct searches to further test the SM, are of great interest. 



This article presents a measurement^ of M\y using data corresponding to a total integrated 
luminosity of 4.3 fb _1 collected from 2006 to 2009 with the DO detector^ at the Fermilab 
Tevatron pp collider. This measurement uses W — > ev events, with electrons in the central 
calorimeter (CC, with \q\ < 1.05, where r/ is pseudorapidity) . The liquid-argon calorimeter 
of the DO detector provides stable and accurate electron energy measurements. The energy 
resolution for an electron in CC at 45 GeV is 4.2% on average □. 

The data set used in this analysis has been recorded with increased instantaneous luminosity 
(almost 3 times higher than in a previous DO analysis^. Higher instantaneous luminosity 
increases the number of additional pp interactions (pile-up) that contaminate the detector. They 
bias the electron energy response and complicate the modeling of the electron reconstruction 
efficiency. Therefore, new developments of the analysis techniques are necessary. 

2 Analysis strategy, event reconstruction 

We reconstruct two vector variables in the plane transverse to the beam direction from a W — > ev 
event, namely, the electron transverse momentum (pt) and the transverse momentum of the 
hadronic recoil (ut) that balances the transverse momentum of the W boson. 

The electron energy is reconstructed as a sum of the energies of calorimeter cells inside the 
electron reconstruction cone, while the direction of the electron is given by the track in the inner 
detector that matches spatially to the calorimeter cluster. The electron energy measurements 
are corrected for the energy loss due to uninstrumented material in front of the calorimeter. The 
correction is derived using detailed first-principle simulation. The material budget is determined 
from a fit to the longitudinal energy profile in the electromagnetic (EM) calorimeter. The gains 
of the readout cells of the EM calorimeter are calibrated using Z — > ee events taking the world 
average Z boson mass^ (Mz) as reference. 

The ut is reconstructed by a vectorial sum of the transverse energies of all the calorimeter 
cells outside the electron reconstruction cone. The longitudinal component of the hadronic recoil 
cannot be determined due to the limited pseudorapidity coverage (|?y| < 4.2) of the calorimeter. 
Therefore, the neutrino longitudinal momentum, which is required to reconstruct the invariant 
mass of the W boson, cannot be determined. 

From and ut, we can calculate three transverse observables for M\y extraction: the 



transverse mass of the W boson (mj- = y (p T +Ifi T ) 2 — (prp +^£ T ) 2 ), the electron transverse 
momentum (p T = \pt\), and the missing transverse energy ($ T = \ — p^ — ut\) due to the 
neutrino transverse momentum. 

A fast Monte-Carlo (MC) model is developed to generate a series of templates for the above 
three observables based on different Mw hypotheses. The Myy is determined, specially for each 
observables, using a binned likelihood fit of the predicted templates to the data. 

For the W —> ev event selection, we require an electron in CC with p T > 25 GeV. The event 
is required to satisfy lfi T > 25 GeV, ut < 15 GeV, and 50 < vtlt < 200 GeV. The requirement 
on ut is made to constrain the transverse boost of the W boson, since the transverse boost 
of the W boson degrades the sharpness of the Jacobian edge in the p T distribution. However, 
this treatment also translates certain uncertainties from the hadronic recoil modeling to the p T . 
There are 1 677 394 candidate W — > ev events after selection. 

The Z — > ee events are the control sample for tuning the fast MC, such as the electron 
energy scale and the hadronic recoil model. The Z — > ee events are selected by requiring two 
electrons both with p T > 25 GeV. Events are also required to have ut < 15 GeV to constrain 
the transverse boost of the Z boson, and 70 < m ee < 110 GeV, where m ee is the invariant mass 




"The electron energy resolution depends on the incident angle (or rf) which characterizes the material budget 
in front of the calorimeter. 



of the electron pair. There are 54 512 Z — > ee candidate events with both electrons in CC, which 
are used for most of the model tuning. Events allowing one electron in the end calorimeter (EC, 
with 1.5 < < 2.5) are only used for measurements of the electron reconstruction efficiency. 

3 Fast Monte-Carlo model 

The fast Monte-Carlo (MC) model for template generation has to simulate W and Z boson 
production and decay, the electron energy response, the hadronic recoil, the underlying events 
contamination, the electron reconstruction efficiency, and the background. 

3.1 Boson production and decay 

The boson production and decay are simulated using RESBOS ^ event generator combined with 
PHOTOS^. RESBDS is a next-to-leading order event generator including next-to-next-to-leading 
order logarithm resummation of soft gluons. PHOTOS generates up to two final state radiation 
(FSR) photons. Parton distribution functions are described using CTEQ6.6HD The boson 
transverse momentum prediction in RESBOS is dominantly determined by the nonperturbative 
parameter^ g2 . The g 2 value^ 0.68 ± 0.02 GeV 2 is used. 

3.2 Electron energy response 

The electron energy response is modeled by firstly modeling the energy responses that are not 
a linear function of the electron true energy. Then, we assume the rest of the energy response 
is a linear function of the electron true energy, fit to the Z — > ee data sample to determine the 
scale. 

The energy loss correction, as one of the non-linear energy responses, is applied to the data. 
There are also certain non-linear energy responses due to the high instantaneous luminosity to 
be modeled in the fast MC. 

One of them is the reduction of EM calorimeter gain due to a high voltage (HV) drop 
caused by a large instantaneous pile-up energy deposition. A large current that flows through 
the resistive coat of the HV pads of calorimeter cells creates a reduction of the HV. The HV 
supplies are connected at both ends of the CC modules (at \n\ = 1.2). Therefore, HV drop 
is larger for cells at small |r/| than for cells at large \rj\. This gain loss is modeled as function 
of instantaneous luminosity and detector r] in the fast MC. The EM calorimeter calibration 
at the cell level applied to the data is done in the absence of the knowledge of the HV drop. 
Certain imperfections in such a calibration are expected. Thus, an addition model of the residual 
miscalibration as a function of detector r] is introduced in the fast MC. 

Another non-linear energy response is understood as an effect of the electron reconstruction 
cone. The electron energy is reconstructed as a sum of energies deposited in a cone consisting of 
13 calorimeter towers. Not only the electron deposits its energy into this cone, but also some of 
the hadronic recoil, pile-up, and spectator parton interactions. The energy deposition from the 
latter sources does not come from the true electron but is reconstructed as part of the electron 
energy in the data. This additional energy contribution is modeled in the fast MC as a function 
of instantaneous luminosity, r], u\\ (the ut projection to the electron direction), and SET (the 
scalar sum of the transverse energy deposited all over the calorimeter with cells in the electron 
reconstruction cone excluded). 

After modeling of the non-linear responses, the linear response is modeled as E = a ■ 
(Etrue — 43 GeV) + (5 + 43 GeV, where, a is the energy scale, f3 is the energy offset, and 43 GeV 
is an arbitrary offset introduced technically to improve the stability of the fit for a and f3. 
The parameters a and /3 are determined by a template fit to the m ee versus fz distribution 
of the Z — > ee events, where fz = {E\ + E 2 ) • (1 — cos uj)/m ee , E\ and E 2 are energies of 



the two electrons, and u is the opening angle between the two electrons. The a and /3 are 
determined separately for four instantaneous luminosity sub-samples, and are consistent with 
each other, as shown in Figure [1] (c). After the electron energy scale tuning, a Z boson mass 
fit returns Mz = 91.193 ± 0.017(stat) GeV, which is in good agreement with the world average 
(M z = 91.188 GeV). The M z fit is shown in Figure □ (a) . 




Figure 1: (a) The dielectron invariant mass distribution in Z — > ee data and from the fast MC, (b) the \ plot of (a), 
and (c) the fitted scale and offset 1-sigma contours in bins of instantaneous luminosity (in units of 10 32 cm _1 s _1 ). 



3.3 Hadronic recoil 

The hadronic recoil (ut) reconstructed from the data contains contributions from the hadrons 
recoiling against the W boson, and additional components that are (mostly) independent of the 
W boson boost. These additional components include spectator parton interactions, pile-up, a 
small part of all the components above that enter the electron reconstruction cone, and the FSR 
photons outside the electron reconstruction cone. 

All these components are modeled separately in the fast MC, and a vectorial sum of them 
gives the ut- Free parameters are reserved to tune the recoil model. The hadronic response 
and resolution are tuned using the mean and width of the r/; m b distributions, respectively, of 
the Z — > ee events in bins of , where ?7i m b denotes the projection of the sum of the dielectron 
transverse momentum and ut vectors on the axis bisecting the dielectron directions in the 
transverse planed 



3.4 Efficiency 

The kinematic dependence of the electron reconstruction efficiency sculpts the distributions of 
the observables used to measure the My/. The pile- up due to high instantaneous luminosity 
and the hadronic recoil are the two major sources that contaminate the electron reconstruction 
window and give inefficiency in the electron reconstruction. The effects of pile-up can indirectly 
introduce kinematic dependence of the efficiency. For instance, under a given pile-up contam- 
ination, a high energy electron can be more easily identified than a low energy electron. The 
hadronic recoil contamination depends on the amount of hadronic activity and also the relative 
orientation of the hadronic recoil with respect to the electron. 

The kinematic dependence of the efficiency is determined and modeled in the fast MC in 
two steps. In the first step, we extract the efficiency dependence from a high statistics detailed 
GEANT^^ MC simulation (full MC) of the W — > ev and Z — >• ee events generated using the 
PYTHIA^l^ event generator, and model it in the fast MC. The full MC is overlaid at the cell level 
with a dedicated pile-up collider data sample which is weighted according to the instantaneous 
luminosity distribution (with random bunch crossings) of the data set used in this analysis. 



This step simplifies the modeling of the complicated correlations among the various efficiency 
dependencies with a data-based determination of the pile-up impacts. The second step is to 
extract the dependence of the efficiency on the major variables (pj., Uu, SET, instantaneous 
luminosity, etc.) from the data, and compare them with those from the full MC. Excellent 
agreement is found between the full MC and collider data. 

3. 5 Backgrounds 

Backgrounds in the W boson candidate sample modify the shapes of the distributions of the 
three observables. The major backgrounds are Z — > ee events where one of the electrons 
escapes detection, multijet events where a jet is misidentified as an electron with lfi T arising 
from misreconstruction, and W — > tv — > evvv events. The fractions of the backgrounds in the 
W boson candidate sample are 1.08% for Z — > ee, 1.02% for multijet events, and 1.67% for 
W — > tv — > evvv events. The impact of the uncertainties in the background model on the My/ 
measurement is found to be small. 

4 Results and outlook 

4-1 Systematic uncertainties 

The systematic uncertainties of the My/ measurements are listed in Table HJ They are divided 
into two categories: from experimental sources and from boson production and decay modeling. 

Among the experimental aspects, the uncertainties from electron energy calibration, electron 
energy resolution model, and hadronic recoil model are driven by the limited statistics of the 
Z — > ee control sample. The shower modeling systematic uncertainties reflect the uncertainties 
in the amount of uninstrumented material, and the energy loss systematic uncertainties arise 
from the finite precision of the simulation of electron showers based on a detailed model of 
the detector geometry. The systematic uncertainties of electron calibration, electron resolution, 
electron reconstruction efficiency, hadronic recoil model and backgrounds are determined by 
varying the corresponding parameters within the statistical uncertainties of their measurements. 

The uncertainties due to boson production and decay modeling are dominantly due to the 
PDFs. In principle, the transverse observables used in the My/ measurement are insensitive to 
the uncertainties of the (longitudinal) PDFs. However, our requirements on the lepton pseudora- 
pidity (|r/| < 1.05) is not invariant under longitudinal boosts. Changes in the PDFs can modify 
the shapes of the transverse observables in the presence of the pseudorapidity requirements. 
The PDF uncertainties are prop agated to My/ by generating ensembles of W boson events 
using PYTHIA with CTEQ6.1^^I The QED uncertainties are estimated by comparing PHOTOS 
to WGRAD 1 - and ZGRAdQ^ event generators, which provide a more complete treatment of elec- 
troweak corrections at the one radiated photon level. The uncertainties from boson transverse 
momentum modeling is determined by propagation of the uncertainty of the g2 parameter. 

4-2 Results 

The value of My/ is extracted by fitting templates of the three observables (m^, pj,, and^ T ) 
generated by the fast MC to the distributions from the collider data. The fitting results are 
shown in Table [2j together with optimized fit ranges, for the three observables. Figure [2] shows 
the distributions of the three observables in data and the comparison with templates from fast 
MC for the best fit My/. During the tuning of the fast MC to describe the collider data, an 
unknown constant offset is added to the My/ values returned from the fits, which is the same for 
the three observables. This enables the tuning of the fast MC to be done without the knowledge 
of the final results. 



Table 1: Systematic uncertainties of the Mw measurement. 





AM W (MeV) 


Source 


rriT 


Pt 




Electron energy calibration 


16 


17 


16 


Electron resolution model 


2 


2 


3 


Electron shower modeling 


4 


6 


7 


Electron energy loss model 


4 


4 


4 


Hadronic recoil model 


5 


6 


14 


Electron efficiencies 


1 


3 


5 


iJfiCK^i OUI1U.S 


2 


2 


2 


Experimental subtotal 


18 


20 


24 


PDF 


11 


11 


14 


QED 


7 


7 


9 


Boson pt 


2 


5 


2 


Production subtotal 


13 


14 


17 


Total 


22 


24 


29 



Table 2: Results from the fits to data. The quoted uncertainty is solely due to the statistics of the W boson 

sample. 



Variable 


Fit Range (GeV) 


M w (GeV) 


xVdof 


rriT 


65 < m T < 90 


80.371 ±0.013 


37.4/49 


Pt 


32 < p e T < 48 


80.343 ± 0.014 


26.7/31 




32 < $ T < 48 


80.355 ±0.015 


29.4/31 



Combining the results from txit and p T methods using the BLUE^W method, we obtain the 
final result of the 4.3 fb _1 measurement: 

M w = 80.367 ± 0.013 (stat.) ± 0.022 (syst.) GeV 
= 80.367 ± 0.026 GeV. 

The result from the^ T method is not used in the combination. In the combination, we assume 
100% correlation for those uncertainties that are nonstatistical in nature, such as the QED 
uncertainties, to protect them from being decreased. However, with this protection, the BLUE 
combination gives a sizable negative weight for the Mw value from the^ T method which has 
relatively larger uncertainties. The interpretation is that the central values of M\y from the 
other more precise methods (m-r and p T ) have also fluctuated apart from the true value of 
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Figure 2: The (a) mr, (b) p^, and (c)^ T distributions for data and fast MC simulation with backgrounds, together 
with the x plots comparing data and fast MC. The fit ranges are indicated by the double-ended horizontal arrows. 



Myy in the same direction as the less precisely method. Given that it is our protection that 
introduces the negative weight, and the contribution to the combined precision from the Tfi T 
method is negligible, we decide to only use the my and methods in the above combination 
to avoid the potential bias. 

We combine our result with the earlier DO measurement^ to obtain the DO 5.3 fb _1 result: 

M w = 80.375 ± 0.011 (stat.) ± 0.020 (syst.) GeV 
= 80.375 ± 0.023 GeV. 

The precision achieved is the same as the previous world average. 

The combination with all previous measurements and the recent CDF measurement^^ gives 
the new world average^S 

M w = 80.385 ± 0.015 GeV. 

The indirect constraints^ on the Higgs boson mass based on this new world average is M# = 
94^24 GeV, with an upper limit of 152 GeV at 95% C.L.. A summary!^ of the measurements 
of the M\y and their average is shown on Figure [3] (left). The world averages of Mw and Mt 
measurements with still allowed regions after direct searches of the Higgs boson are shown in 
Figure [3] (rig ht)Hl. 



Mass of the W Boson 

Measurement M„ [MeV] 

CDF-0/1 ■ — ' • ■ 80432 ±79 

D0-I ' • ■ 80478 ±83 

D0-II ■ — • — ■ 80402±43 

CDF-Ilaart 80387±19 

D0-II (4.3 b _l ) 80369±26 

Tevatron Run-0/l/ll 80387 ±16 

LEP-2 — 80376 ±33 

World Average 80385 ± 1 5 
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Figure 3: (left) Summary of the measurements of the W boson mass and their average as of March 2012. (right) 
1-sigma contour in the (Mt, Mw) plane representing the current world averages of Mt and Mw measurements. 
The green bands show the possible Higgs boson masses that are not excluded by direct searches. 



4.3 Outlook 

DO has another ~ 5 fb^ 1 collider data to be analyzed. Including this last data set and still only 
using electrons in CC, if we assume the uncertainties from boson production and decay model 
would be unchanged, the precision of M\y from the total ~ 10 fb _1 DO data is expected to be 
~ 18 MeV. 

If the PDF uncertainties would be reduced by a factor of two by using new PDF sets (e.g. 
Ref.^31) including constraints from W charge asymmetry measurements, and also using electrons 
in EC to extend the pseudorapidity coverage, the above precision can be improved to ~ 16 MeV. 
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